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^ ! Abstract 

We have reviewed recent developments of the theory of the impact for macroscopic elastic ma- 
terials. This review includes (i) standard theories for the normal impact and the oblique impact, 
(ii) some typical approaches to simulate impact problems, (iii) and an example of our simulation 
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to clarify the mechanism of anomalous restitution coefficient in an oblique impact in which the 
restitution coefficient exceeds unity. 
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I. INTRODUCTION 



Impacts are common in nature. Besides microscopic impacts for atoms and molecules, 
there are plenty of examples of impacts for macroscopic materials. To control impacts is 
important in ball games in sports and many processes in industrial plants. The impacts 
of such the macroscopic materials cause complicated processes and eventually they become 
inelastic. 

Studies of inelastic impacts are aimed to clarify the relation between the pre-collisional 
state and the post-collisional state. Since there are huge number of degrees of freedom in 
macroscopic materials, it is difficult to follow all the processes of energy transfers induced 
by the impact. In this situation, we need to introduce simple quantities to characterize the 
inelastic impact of macroscopic materials. 

For this purpose, the coerhc.ent of restitution e is w,de, y u S edQ Q Q Q Q, which is 
defined by 

v c • n = -ev c ■ n, (1) 

where v c and v c are respectively the relative velocity at the contact point of two colliding 
materials before and after the collision, and n is the normal unit vector of the tangential 
plane of them (Fig.l). Although many text books of elementary physics state that e is 
a material constant less than unity, it has been confirmed that e decreases as the impact 
veloc.ty _|MflQ. For example, the dependence of e for the low impact velocity 
are theoretically treated by the quasi-static theory 0, 0, [lfl, U, 12]. 

Rebound processes depend on the impact angle. Therefore, we also introduce the coeffi- 
cient of tangential restitution j3 as 

v c ' ■ t = -/3v c • t, (2) 

where t is the unit tangential vector (Fig.l). (3 is a function of the incident angle 7 which 
is defined as 7 = axctax\.{yt/v n ) with y n = v c • n and v± = v c • t. It is believed that possible 



values of (3 lie between -1 and 1 
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. The relation between the impact 



speed v of the center of mass and v c at the contact point is given by 

v c = v — Ruj x n (3) 

for hard spheres, where R and uj are the radius of the sphere and the angular velocity 
of the sphere, respectively. The phenomenological theories of the oblique impact have been 




FIG. 1: A schematic picture of collision between a sphere and a flat wall. 



The impact process is important in granular physics. Although we believe that the 
static interaction among contact grains can be described by the contact mechanics of elastic 
materials, we little know the dynamical part of contact mechanics. Thus, the distinct 
element method which is the most popular method for the simulation of grains contains 



many unknown parameters. 



the basis of granular physics. [2c 



herefore, the theory of impact for elastic materials gives 
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While e has been believed to be less than 1 in most situations, we have recently recognized 



that e can exceed 1 in oblique impacts 



27L l28j |. In particular, Louge and Adams 



reported that e increases as a linear function of the magnitude of tan 7 in the oblique impact 
of a hard aluminum oxide sphere on a thick plate with the incident angle 7. In this case, 
Young's modulus of the wall is 100 times smaller than that of the sphere in the experiment. 
Thus, the physics of impact is one of interesting subjects in current statistical mechanics. 

The organization of this paper is as follows. In the next section, we explain the current 
understanding of normal impacts including the quasi-static theory, the effect of radiation of 
sounds, and the effect of the plastic deformation. In section III, we summarize the standard 
treatments for oblique impacts which include Walton's argument 14 1 and the theory by 
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Maw et a/.Q| In section IV, we will introduce three typical approaches for simulation of 
impacts of elastic materials. In section V, we briefly explain the recent analysis to explain 
the anomalous behavior of the restitution coefficient exceeds unity. Section VI is the short 
summary of this paper. 

II. THE CURRENT UNDERSTANDING ON NORMAL IMPACT 

This section is devoted to summarize the current understanding on a normal head-on 
collision of e 



astic spheres or a collision between a sphere and a wall. From the quasi- 

B 
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121 ] . we believe that the restitution coefficient e decays as 1 — 
e oc v 1 ^ 5 with v = v n for small impact velocity. The agreement between the theory and 
experimental results is fair, j^, 0, ^ For high speed impacts, the plastic deformation of elastic 
particles is dominant mechanism to determine the post-collisional processes. Experimental 
results support the theoretical prediction e cx v -1 / 4 .)^! 2^ At present, we do not have 
any appropriate theory for the finite impact speed but below the threshold at which plastic 
deformation takes place. 



A. Mechanism of inelastic collision 

Inelasticity arises from the transfer of translational kinetic energy to internal degrees 
of freedom. The dominant dissipative processes for low speed impacts are viscous effects 
of elastic materials, the heat conduction and the sound emission into or out of the elastic 
materials. When a local deformation with finite speed takes place in a macroscopic material, 
the system is excited to a nonequilibrium state and after that it is relaxed to an equilibrium 
state. It is not easy to specify the microscopic origin of viscous effects or the relaxation 
process, because such systems have numerous number of degree of freedom. However, at 
least, inelastic scatterings of phonons and excitation-radiation processes in electronic states 
are two major sources of the dissipations. 

It should be noted that the coefficient of restitution of one- dimensional rods is insensitive 
to the impact speed but depends on the ratio of lengths for two colliding rods. 
On the other hand, the restitution coefficient strongly depends on the impact speed and 
Poisson's ratio in higher dimensional impacts. For example, the two-dimensional simulation 
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by Gerl and Zippeliusj33j shows that inelasticity increases as Poisson's ratio increases. 



B. Outline of theory of a quasi-static impact 

In this subsection, we present the outline of quasi-static theory for a normal impact of 
elastic spheres. 

The theor, * based on the contact t he 0ry of elastic sphe.es devetoped by He rt ,Q 
Hertzian contact theory predicts that the radius of contact a and the elastic compress force 
F e i for two spheres with radii R and R! are respectively given by 

rW„ RR ' r, ( RR' 

a = Ftr D : F d = (4) 

I R + R' J D \R + R'J K J 



where h is the length of compression, and D = + 



3ri-v i i -u' 2 \ 
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with Young's modulus E, E' 
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38] 



and Poisson's ratio u, v' for two contact materials. 

Let us consider a low speed impact of two spheres. In the limit of low speed, the nonequi- 
librium processes may be suppressed, and the collision can be treated as an elastic process. 
The energy conservation can be read 
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dh Y ,5/2 2 4 RR' 
m ^[dt) +Kh =me ^' k= 5D^rTr^ (5) 

where m e ff = mm' /{m + m') is the reduced mass of two spheres with masses m and m' . 
The maximum deformation is easily obtained as h = (m e ///K) 2 ' 5 u 4 ' 5 . The contact time t c 
of the collision which is two times as large as the time needed to reach h is given by 

ic = 2 (!%) 1/ 7 1 ^ f = C( (!%) 1/5 (6 ) 

\ K Z V j JO VI- X 5 / 2 V K V J 

where c t = 4 v /7rr(2/5)/5r(9/10) ~ 2.94 with the Gamma function T(x). 

The above treatment predicts e = 1, because the process does not contain any dissipation. 
Kuwabara and Konop] assume the existence of Rayleigh's dissipation function for elastic 
solids, and write down the dissipative force as 



/ RR' Vhdh 
Ff = ~lRTK^H (7) 

for two contacted spheres. Here, D corresponds to D for elastic contact, which can be 
represented by viscous parameters of Rayleigh's dissipation function. The magnitude of the 



viscous parameter D can be measured from an experiment of sound attenuation. Adding 
this viscous term, the dimensionless form of equation of motion becomes 

5 5^ / v X 1 / 5 

x + r)^x + -a: 3 / 2 = 0; 77 = -k [-^tj-^ (8) 



where k = (A/hD)jRR'/( K R-\-R l ). Here, we nondimensionalize the variables in terms of 
h = h x, t = (h /v)T, and thus x = dx/dr. Thus, the problem is reduced to obtaining 
e = dx/dr at r c = t c (v/ho) under the initial conditions x = and x — 1. From (JHJ), it is 
easy to obtain 

e 2 - 1 = -2r) r Jxx 2 dr. (9) 
J 

Since the exact evaluation of the integral is impossible, we may replace x in the above 
equation by the solution of elastic equation. Using this approximation we obtain 

5 / \ 1/5 

e ~ 1 - 1.00977 = 1 - 1.009 x -R ( — - , (10) 

2 \K 3 m eff 2 J 



where the numerical constant comes from (4/5)5(3/5,3/2) ~ 1.009 with the beta function 
B(x,y). Thus, the coefficient of restitution for the low speed impact is believed to decay 
1 — e oc v 1 ^ . The result can be obtained in different contexts. P. u\ 

Here, we briefly summarize the two-dimensional counterpart of normal impacts. For 
impacts of an elastic disk on a rigid wall, we do not have reliable argument. The total force 
of elastic force and dissipative force may be given by 

tot ~ \n(AR/h) T0 \n(4R/hY [ } 

where = E/(l — u 2 ) and tq represents the time scale for the dissipation of the small 
deformation. Replacing the logarithmic term as a constant correction, the equation for 
elastic motion can be solved. Thus, we may evaluate the contact time t c as 



t c ~— Win— 12 
c V v 

where c = \J E*/ p and p are the compressive sound velocity and the density, respectively. 0, 
39| From the comparison of eq. ([T2|) with a two-dimensional simulation, we have confirmed 
that the above estimation is quantitatively correct. 39] Here, we adopt a bold approximation: 
h max ~ vRJ \pj 'E*)\n(4R/h) ~ vR/c. Including the dissipative force this approximation 
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gives 

1 27tE * (i^ 
e~l . (13) 

pRc^ln(4c/v) 

The correctness of this expression has not been confirmed. 
C. The effect of radiation of elastic waves 

So far, we do not have any theoretical argument to estimate the restitution coefficient 

nr 

as a result of radiation of elastic waves. On the other hand, Miller and Pursev|40l. 141 1 
calculated the averaged power radiated per unit area by the normal oscillating contact of 
circular region on a semi-infinite isotropic elastic material. Although the situation is a little 
different each other, we may apply their calculation to estimate the energy loss by radiation 
of elastic waves. 

Miller and Purseyj^, l^j] obtained the powers radiated in the compressible (W c ) and 
shear waves (W s ) as 



w c = — 1 ^ y o de{Q x {e)Y sine (14) 

W s = 1 p / d6\e 2 (6)\ 2 sm6, (15) 



4p Jo 

where 6 is the polar angle from the axis of symmetry, Pq = F e i/(iia 2 ) is the average com 
pressive force, and 



F o (sin0) Fo(Qsm.v) 



F (x) = (2x 2 - C 2 ) 2 - 4xV(^ 2 " C 2 )(^ 2 - 1), (17) 



7T 



+ !/)(!- 21/) vr 2p(l + i/) 

^ ^(l-^) ' k2 ~d^E— (18) 

Here we omit the power radiated in terms of surface waves which plays an important role 
in the paper by Miller and Pursevfiol. [il l, because they can be absorbed as the numerical 
factor and the treatment for spherical surfaces is not obvious. The integrations of W c and 
W s are possible once we specify the material. For example, in the case of the material with 
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v — 1/4, i.e. £ = a/3, the total amount of power is given by 

W — W c + W s — 1.579 



a 4 fc! 2 t c C 4 P 2 



4p 



The energy loss during collision is thus given by 

rtc/2 _ 5.337 



E 



loss 



(20) 



(21) 



Similarly the argument in the previous section, the integral in the right hand side can be 
evaluated as J * c/2 ' dtF% = 5 v ^Fr(8/5)/(8r(21/10))/? 2 / 5 m^w 11 / 5 in the elastic limit. Thus, 
we obtain 

k 6 / 5 ™ 4/ 5 .... 

(22) 



K 6/5 m 4/5 

E lnss ~ 0.5827 °/V 3/5 . 



pE*/ 2 



From the relation E, 



loss 



[l-e 2 )m eff v 2 /2 



0.5827- 



m e ffV (1— e), we finally reach the new relation: 

„6/5 



1/5 



(23) 



pm x J f )EW 

for i/ = 1/4. The discussion here is the rough evaluation of the restitution coefficient based 
on the energy loss by radiation of elastic waves. Although the numerical constant in eq. (|23|) 
is meaningless, we expect that the scaling relation can be used in realistic situations. It is 



interesting that the restitution coefficient e o 
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42 1 derived 1 — e oc v 3//5 from the theory of Miller and 



as that of the quasi-static theory. [6 
It should be noted that Hunter 
Pursey. j^j, 41] The difference comes from the followings. Hunter 
dissipation is obtained from 



j eys 1 — e oc v 1 ^ which is the essentially same 
13 



42j assumes that the energy 



W 



H 



tc dt™ 2 F d {t)^p-. 
at 



(24) 



where u(t) is the mean surface displacement. Although Hunter 42j adopts the result of Miller 
and PurseyQ, |^| for u, this is only the fraction of the radiation. The above treatment 
presented here is more reliable. 

However, the analysis presented in this subsection is still prematured. Insufficient parts 



of the analysis are as follows: (i) The theory by Miller and Purseyj^j, |4l| assume the 



constant pressure in the contact area, but Hertzian contact theory predicts the distribution 
of pressure, (ii) Miller and Purseyj40j, |4j| discussed the radiation of elastic wave for semi- 
infinitely large region, but the actual contacted spheres have curvature and finite volume. 
Therefore, we may need more systematic treatment for the radiation of elastic waves. 
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D. Impact with plastic deformation 



When the impact speed is large enough, the elastic description is no longer valid but we 
have to consider the effects of plastic deformation. The restitution coefficient drastically 
decreases when the plastic deformation occurs. Following the argument by Johnson ;30], we 
review the argument of the restitution coefficient for impacts with plastic deformations. In 
the argument in this subsection, we neglect numerical factors. So the result is basically valid 
for collisions for two spheres with equal radius R and the density. 

Let us evaluate the yield pressure. From Hertzian contact theory the maximum pressure 
of compression exists at the center of contact and its expression is given by 

2\ 1/3 

2vra 2 V tt 3 R 2 ~ 



3F e i (QFeiE* . 
Po = ^7 = _ 3D o • (25) 



Thus, the compressive elastic force at the yield pressure p = (po)y becomes 

Fei = ^(Pofv, (26) 

where (po)y is about 1.6 E at Mises' condition. [30!] 

Hertzian contact theory gives F e i ~ E*R 1 / 2 h 3 >' 2 , while the maximum deformation is 

; 2 \2/5 



h * ~ (27) 



from mv 2 ~ F e [h ~ E^R 1 ^ 2 ^^ 2 , where m is the mass of each sphere. Substituting this h* into 
Hertzian contact theory with setting F e i = Fy, we balance the force with Fy ~ ^-(po)y- 
Then we obtain 

/ F \ 4 / 5 



(Po)y - [j^) {mv 2 y\ (28) 

which is the condition for yield stress. 

Let us discuss the effect of plastic deformation to the restitution coefficient. From 
Hertzian contact theory there is a relation between the compression and the radius of 
the contact h ~ a 2 / R. Therefore the deformation for the plastic deformation becomes 
/i* ~ F*/(E*a*) because of F* ~ E*y/ Rhh, where a* and are respectively the radius 
of contact area and the compressive force for the plastic deformation. Thus, the internal 
energy stored during the deformation may be evaluated as W ~ F^h^ ~ F 2 /(a*E*). On 
the other hand, we have F* = na^p^ with the contact pressure pd- Since the stored energy 
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is released as the kinetic energy in impact processes, the kinetic energy for rebound is given 
by K r ~ W> ~ affi/E.. 

On the other hand, we can estimate the work needed for compression as W = (l/2)mv 2 ~ 
Jq* dhF e i ~ Jq" daF^a/R. Assuming that F e \ = ira 2 pd is kept during the impact, we obtain 

u/ f a * irp d a 3 da p d a* 4 2 

W = / ~ ~ mv . 29 

Jo R R v ' 



Thus, we reach 

2 v' 2 p d R 

Equation ()30|) can be rewritten as 



(30) 



2 (p d R 3 \ l/4 

e ^^l^J (31) 

and thus, we finally obtain 

e oc v- 1 ' 4 . (32) 
This expression recovers experimental results. 

III. CURRENT UNDERSTANDING OF OBLIQUE IMPACT 
A. Walton's argument 

To characterize the oblique collision, Walton introduced three parameters: Q, Q the 
coefficient of normal restitution e, the coefficient of Coulomb's friction /x, and the maximum 
value of the coefficient of tangential restitution (3 . The expression is given by 

^|-l-Ml + e)cot 7 (l + Hf) ( 7 > 7o ) 
[Po (7<7o), 

an gle, and m, R, and / are mass, radius and moment of inertia 
13!. llJ| Experiments have supported that his characterization ade- 



where 70 is the critica 
of spheres respectively. 



quatelv capture the essence of binary collision of spheres or collision of a sphere on a flat 
plateQ, 16, 17, 3- 



The derivation of the first equation of Walton's expression is simple. When there is a 
slip, the friction coefficient satisfies 

|n x J| = /i(n ■ J), (34) 
10 



where J = m(v' — v) is the impulse. Let us write J as the following form: 

J = A(n ■ v c )n + B[v c - (n ■ v c )n], 



(35) 



where A and B are constants to be determined. From the definition of e and (3 with the aid 
of I(uj' — u>) = — R(n x J) we write 

m(v' c -v c ) = A(n-v c )n+(l + ^-)5[v c -(n-v c )n]. (36) 

Thus, from the projection to the normal direction we obtain A = — m(l + e) and B = 
(1 + mR 2 /I)m(l + (3). On the other hand, through the relation |n x J| = B\n x v c | = 
yu(n • J) = fj,A(n ■ v c ), we can rewrite J as 

J = — m(l + e)(n • v c )n + /im(l + e) cot 7[v c — (n • v c )n]. (37) 

Thus, we obtain the first expression of eq.(J33J). 

1.5, , , , ^-^ 




FIG. 2: Comparison between Walton's argument and our numerical result. Here = tan 7 and 
^2 = —vt'/v n . Parameters for the simulation are the same as those in ref. . 

In spite of its simple form and its simple derivation, Walton's expression is useful to 
characterize oblique impacts. However, we do not know how to determine /3q and 70. In 
addition, it cannot be applied to impacts for very small 7 in which (3 cannot be a constant, 
because it should not be contradict with the normal impact at 7 = 0. The discontinuity of 
differentiation of f3 at 7 = 70 is also unnatural. To know the details of physics of the oblique 
impacts we should introduce alternative theory. 
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B. Theory of Maw et al. 

More systematic treatment for oblique impacts is developed by Maw et al. 
complete explanation of their theory is long because of its complicated structure. Here, we 
only summarize the result of the theory. 

They assume that the stiffness of normal compliance for restitution changes from k to 
k/e 2 where k is the normal stiffness for compression. They do not discuss the origin of the 
normal restitution coefficient e which is assumed to be constant and its effect for the stiffness 
explicitly. 

They also indicate that there are three region depending on the incident angle 7. Ac- 
cording to their theory [ljj], (3 for the impact of a circular disk on a wall can be represented 

by[ 



(i) 1/ficr 2 < cot 7 with a = y(2 - u)/(2(l - v))\ 

'Qti 7T 



(3 = — cosuti — fi^-e 

Pz 



l + cos(— + |(l-e" 1 )^) cot 7, (38) 



(ii) p x /p z jj,(l + e) < cot 7 < 1/jia 



2. 



(3 = — COSUj(t 3 — t 2 ) — /J^{cOStu(t3 — t2) — COsVlt 2 COS Uj(t 3 — t 2 ) 

Pz 

H — sin Vtt 2 sin oj(t s — t 2 ) + e + cos fi£ 3 ] cot 7, (39) 

UJ 



(iii) cot 7 < p x /p z fx(l + e): 



= -l + /i§^(l + e)cot 7 , (40) 

Pz 



where /1 is the coefficient of friction, /3 X = 3 and (3 Z = 1 are constants for the impact of the 



disk on the infinite wall. Vt and uj are respectively ir/2to and (ir/2ato)y (3 X / (3 Z} where to is 
the time for compression. The time £1 is the transition time from initial stick motion to slip 
motion which is determined by 



1 v x (0) Q sin ujt 



- = 1 < ti < t 



I^U^OI _ I o" 2 fiv z (0) uj sin Qti , , p 

^F z {ti) ~ 1 t;,(0) fl sina^ _ ^ 

1 aV^W^sin^ + Kl-e- 1 )) ° " * < /j 
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where v x (t) and v z (t) are the tangential and the normal relative velocities at the contact 
point, respectively. While t 2 can be determined by 

'v x (0)/fiv z (0)-(3 x /(3 z \ 



— arccos 

Qt 2 

e 



a- 



7T , 

- 1-e 

2 y 



-i-> 



arccos 



'v x (Q)/fMv z (0)-p x //3 z 



r V X (0) (3 X 

for — — < a— 

v z (0) -*p t 



e Wx(0) As 

tor — — > u— 

v*(0) 



which is the time to start sticking. The time £3 determined by solving eqs. 
is the transition time from stick motion to slip motion: 



Qu x (t 2 ) , , £lv x (t 2 ) . 

COS UJ{ts — t 2 ) 7-r sin UJ{t 3 



a 2 sin 



(42) 
numerically 

(43) 



where u x (t2) is the tangential deformation at time t 2 . By calculating (3 at each value of cot 7 
and interpolating them with cubic spline interpolation method, we can draw the theoretical 
curve. Figure El shows comparison of the theory with our numerical simulation, in which the 

1.5 



0.5 



-0.5. 




0.5 



1 



1.5 



FIG. 3: Comparison between Maw's theory and our numerical result. The situation is the same as 
that in ref.|2fl|. 

agreement is acceptable. 

As is shown, the theory by Maw et al. P, ^| captures the physics of impact processes. 
For large 7, i.e. region (iii), the disk slips on the wall without any rotation or sticking. For 
intermediate 7, that is, the region (ii), the disk slips at first and stick at t = t 2 and slips 
again at t — £3. For small 7 in the region (i), the disk sticks initially and begins to slip 



at t — t\. In this sense, their theory improves some defects of Walton's argument. 
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However, their theory includes some other defects: The final expression for (3 is complicated 

13 



and is required for numerical calculation. It still contains undetermined parameters e and 



\i which are assumed to be constants. As is shown, e and jjl stron 
speed and the incident angle, 
justify their assumption. 
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33, 



y depend on the impact 



33 . |43j Thus, we cannot 



IV. NUMERICAL MODELING 



In general, the success of theoretical approach is limited for nonlinear problems because 
of difficulties of analysis. On the other hand, numerical approach has become standard 
as computers become popular. The advantage of this approach is obvious, (i) We can 
investigate collision processes under idealistic situations, (ii) It is easy to control situations 
to investigate the properties of impact processes, (iii) It is possible to analyze nonlinear 
problems. However, when we restrict our interest in numerical studies of impacts, we still 
do not have any standard technique and the status of such the studies is prematured. 



45 



One of typical approaches for engineers is to use FEM (Finite Element Method). 



46 



44, 



47( There are some standard packages for simulation of impact processes. For exam- 



ple, Lim and Stronge|44 



46J carried out two-dimensional simulation of a transverse collision 



of a cylinder against an elastoplastic half space based on DYNA2D. A three-dimensional 
FEM also exists. (^tJ All of them reproduce experimental results. Since FEM is originally 
proposed for a solver for static problems, they can recover the Hertzian contact theory for 
the static elastic problem. No viscous term, however, is included within FEM. To obtain 
inelastic impacts, FEM usually introduces elastic-plastic deformation or fully plastic defor- 
mation. Therefore, the simulation based on FEM sometimes predicts e oc v _1//4 without 
quasi-static region. |44] We also note that there are many input parameters for FEM. For 
example, it includes at least two yield stresses for the transition from the elastic region to the 
elastic-plastic region, and the transition from the elastic-plastic region to the fully plastic 
region. It also contains at least two friction coefficient as a constant. However, as will be 
shown, the friction coefficient cannot be regarded as a constant. 

Another approach is proposed by Ger , and ZippeliusflQ which is based on the m ode 
analysis of elastic materials. They assume that an isolated disk is in an eigenstate of isother- 
mal vibration of elastic waves. When the disk contacts a wall, the transitions among eigen- 
states is induced by nonlinear effects of interaction between the disk and the wall. 
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The last approach is based on the equation of motion for mass points connected by 
springs. This approach is intuitive and flexible. For example, the other approaches 

may not be applied to impacts with rough surfaces without essential change of algorithms, 
but this approach can include the roughness at surface easily. The roughness at surface 
p.ays c rU c,a> r0 .es for ob.ioue impact, In fact, the coefficient of tangential restitution f, 
becomes -1 without the roughness. [20] The disadvantage of this approach is that the code 
is not fast, and the result strongly depends on lattice structures. 

The last two approaches are adequate for the high speed impact without plastic deforma- 
tion. At present, there is no three dimensional simulations, and two-dimensional simulations 
for normal impacts may not agree with those for quasi-static theory. One of defects for these 
approaches is that the problem is not relaxed to a static state without introduction of local 
dissipation. For later discussion, we focus on the last approach to simulate the oblique im- 
pact of a disk on a flat wall to clarify the mechanism of anomalous behavior of the restitution 
coefficient e. 



V. SIMULATION OF OBLIQUE IMPACT 




FIG. 4: The elastic disk and wall consisted of random lattice. 
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A. Our model 

Let us introduce our numerical model which is a typical example of the last 

approach in the previous section. The discussion in this section is based on our recent 



paper. 



43] 



Our numerical model consists of an elastic disk and an elastic wall (Fig. @J. The width 
and the height of the wall are 8R and 2R, respectively, where R is the radius of the disk. 
We adopt the fixed boundary condition for the both side ends and the bottom of the wall. 
To make each of them, at first, we place mass points at random in a circle and a rectangle 
with the same density, respectively. For the disk, we place 800 particles at random in a 
circle with the radius R while for the wall, similarly, we place 4000 particles at random in a 
rectangle. 

After that, we connect all mass points with nonlinear springs for each of them using 

spring interaction between connected mass 



points is given by 



V®(x) = -k^x 2 + ht^x 4 , i = d(disk), iy(wall), (44) 



where x is a stretch from the natural length of spring, and fcW an d are the spring 
constants for the disk(i=d) and the wall(i=w). 

Inmost of our simulations, we adopt kffi = l.Oxmoc 2 /RHox the disk while fcM = k^ /100 
for the wall, where m and c are the mass of each mass point and the one-dimensional vel ocity 
of sound, respectively. In this model, the wall is much softer than the disk as in ref. 28]. 
We adopt = fcW x 10~ 3 /i? 2 for each of them. We do not introduce any dissipative 
mechanism in this model. The interaction between the disk and the wall during a collision 
is given by F(Z) = £Vo exp(£/)n s , where £ is 300/ R, Vq is £,moc 2 R/2, I is the distance between 
each surface particle of the disk and the surface spring of the wall, and n s is the normal unit 
vector to the spring. 

In this model, the roughness of the surfaces is important to make the disk rotate after 
collisions j^J. To make roughness, the normal random numbers with its average is zero and 
its standard deviation is S = 3 x 10 _2 i? are used for the surface particles of the disk and 
the wall. All the data in this paper are obtained from the average of 100 samples in random 
numbers. 
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Poisson's ratio v and Young's modulus E of this model can be evaluated from the strains 
of the band of random lattice in vertical and horizontal directions to the applied force. 
We obtain Poisson's ratio v and Young's modulus E as v = (7.50 ± 0.11) x 10~ 2 and 
E = (9.54 ±0.231) x 10 3 m c 2 / R 2 , respectively. Q 

In our simulation, we define the incident angle 7 by the angle between the normal vector 
of the wall and the initial velocity vector of the disk(see Fig. We fix the initial colliding 
speed of the disk as |v(0)| = 0.1c to control the normal and tangential components of the 
initial colliding velocity as v t (0) = |v(0)| sin 7 and f„(0) = |v(0)| 00S7, respectively. We use 
the fourth order symplectic numerical method for the numerical scheme of integration with 
the time step At = 10~ 3 i?/c. 




FIG. 5: Numerical and theoretical results of the relation between v&i and e. 



B. Results 

Figure El is the normal restitution coefficient e against = tan 7 for the impact of the 
hard disk on the soft wall. The cross points are the average and the error bars are the 
standard deviation of 100 samples for each incident angle. This result shows that e increases 
as increases to exceed unity, and has a peak around = 6.0. This behavior is contrast 
to that in the experiment by Louge and Adams 3| • 

Let us clarify the mechanism of our results. Louge and Adams [2^ suggest the anomalous 
behavior of e > 1 can be understood by the local deformation on the surface of the wall 
during an impact. They attribute their results to the rotation of normal unit vector of the 
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wall surface by an angle a and derive the corrected expression for e. Thus, we determine 
the quantity of a at each incident angle from the theory of elasticity and calculate corrected 
e. 



y 




- X 



FIG. 6: The schematic figure of a hard disk sliding on a soft wall, x coordinates of both ends of 
the contact area AB are x = x a and x = x\,. 



Figure El is the schematic figure of a hard disk moving from left to right on a wa 
the length of the contact area is I = |a?& — x a \. From the theory of elasticity (35L 
ratio can be estimated as 



X where 
2fi|, this 



x c-x a 1 l-2v 

= 1 — with = — arctan— -, (45) 

/ 7T a*(2 — 2i/) ' V ; 

where v is Poisson's ratio of the wall and /1 is the coefficient of friction. Here tana = 
f(xb)-f(xa) -yfifa fa e deformed shape of the wall f(x) approximated by a parabolic function 
near x c is reduced to 

2^ I X £• a I / a /-1 \ 

tana = . (46) 

2-29 R K ' 

In eq.dlH}, | be evaluated by the simulation data. From our simulation, the 

maximum value of y c is about 0.17-R at ^1 = 1.0. Assuming the disk is pressed in the 
normal direction, we can estimate the contact about 1.1R which is the maximum 

value. Thus, we adopt its half value, 0.55i?, as \x c — x a \. 

The cross points in Fig|7|is \i calculated from eq.(|34"|) and our simulation data against 
each ^1. Figure [7| shows fi has a peak around = 3.0. Substituting this result to eqs.(|4"5~J) 
and ()46|1. we obtain the relation between and tana. 

The restitution coefficient e can be obtained as a function of tana by regarding the 
impact as that on a tilted surface with the angle a. Skipping the derivation, we can write 
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the result asj 

_ e a + ^tana 

6 ~~ 1 - tff tana ' ^ ' 

where e Q is the restitution coefficient defined through 



-e a (v c ■ n a ), (48) 



where n a is the unit normal vector of the tilted slope connecting A and B in Fig. 6. Here 
and ^2 are respectively given by 

¥? = ^~ tana . (49) 
1 l + ^i tan a 

and 

^ = - 3(1 + e a )fi a (50) 



in the two-dimensional situation 



14j. /x a in eq.(|5U|) is given by 



/x + tana 

= i 7 ■ 51 

1 — //tana 



To draw the solid line in Fig0 at first, we calculate tana and \x by eqs.()46|) and 
respectively for each After that, we calculate \&i and by eqs. ()49|) and (|50|l . and obtain 
e by substituting them into eq.(|4T|) for each Although e a is assumed to be a constant 
in ref.J^, here we evaluate e a from our simulation which has the weak dependence of 7. 
The solid line of Figgis eq. ()47j) . All points are interpolated with cubic spline interpolation 
method to draw the theoretical curve. Such the theoretical description of e is consistent with 
our numerical result. As can be seen, the restitution coefficient e depends on the relation 
between /x and 

It should be noted that the behavior of /x as a function of 7 can be understood by using 
a simple phenomenological argument. The solid curve in FigEJis the fitting curve obtained 
from the theory. If we omit the collapse of the roughness by the impact for large 7 close to 
7r/2, both /x and e increase with increasing 7. This s ugg ests that the discrepancy between 
our result in Fig. 5 and the result by Louge and Adams 28] is originated from the difference of 
impact speed. Namely, our impact speed much larger than the experimental one. Although 
we skip the details of derivation, it may be clear that our argument captures the essence of 
physics of the oblique impact. 
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FIG. 7: Numerical and theoretical results of the relation between v&i and fi. 
C. Discussion 



At first, let us discuss the origin of the relation between e and ^i. As was indicated, 
the local deformation of the wall's surface is essential for e to exceed unity. We have also 
carried out the simulation when = 10 x which means the wall is harder than the 
disk. In this case, e takes almost constant value to exceed unity suddenly around = 4.5. 
This tendency resembles the experiment by Calsamiglia et al^27\. This behavior can be 
understood as that the disk is scattered by hard nails distributed on the surface of the wall. 
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FIG. 8: Numerical results of e for the impact of a soft disk on a hard wall. It should be noted that 
= 4.5 is almost the same as the largest 7 in the experiment of ref . (23] - 



Thus, the wall should be much softer than the disk to get smooth increases of e as 



20 



increasing 7. In addition, it is important to fix the initial kinetic energy of the disk. We 
have confirmed so far that e cannot exceed unity when we change 7 with fixed v n [soj] . 

Second, the initial velocity of the disk and the local deformation of the wall are so large 
that the local dissipation in springs and the gravity have not affected our numerical results. 
In addition, we have carried out the other simulation with a disk of 400 particles and a wall of 
2000 particles to investigate the effect of the model size. Although there is a slight difference 
between the results, the data can be reproduced quite well by our phenomenological theory. 

Third, the local deformation of the wall also affects the relation between u and ^1 . In 
early studies, it has been shown that /1 depends on the impact velocity H Q. In our 
simulation, the magnitude of a has a peak around = 3.0. This behavior is interpreted 
as that the local deformation collapses for large v t . The decrease of [i and the friction force 
cause the decrease of e. 

In final, we adopt the static theory of elasticity to explain our numerical results for the 
discussion here. However, it is important to solve the time-dependent equation of the defor- 
mation of the wall surface to analyze the dynamics of impact phenomena. The dynamical 
analysis is our future task. 

In summary of this section, we have carried out the two-dimensional simulation of the 
oblique impact of an elastic disk on an elastic wall. We have found that the restitution coef- 
ficient e can exceed unity in the oblique impact, which is attributed to the local deformation 
of the wall. The relation between /1 and is also related to the local deformation and can 
be explained by a simple theory. 



VI. SUMMARY 



We review the current understanding on inelastic impacts of elastic materials. We explain 
the standard theories for the normal impact and for oblique impacts. We also introduce some 
typical approaches for numerical simulation for impact problems. Although this problem is 
fundamental and familiar in elementary mechanics, the theoretical treatment is prematured. 
In recent finding of the anomalous behavior of the restitution coefficient exceeds unity is an 
example to have potential to be developed as a subject of physics. Unfortunately, most of 
the existing theories are based on engineering idea and they are not so simple and beautiful. 
In some cases, the assumption of the theory may not be valid. For physicists, in other words, 
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there will be a lot of room for improvement of the theory of impact of elastic materials. 
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